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ABSTRACT 

Recent radio and optical large scale jets' observations suggest a two- 
component jet morphology, consisting of a fast central spine surrounded with 
a boundary layer with a velocity shear. We study radiation of electrons accel- 
erated at such boundary layers as an option for standard approaches involving 
internal shocks in jets. The acceleration process in the boundary layer yields in a 
natural way a two component electron distribution: a power-law continuum with 
a bump at the energy, where energy gains equal radiation losses, followed by a 
cut-off. For such distributions we derive the observed spectra of synchrotron and 
inverse- Compton radiation, including comptonization of synchrotron and CMB 
photons. Under simple assumptions of energy equipartition between the relativis- 
tic particles and the magnetic field, the relativistic jet velocity at large scales and 
a turbulent character of the shear layer, the considered radiation can substan- 
tially contribute to the jet radiative output. In the considered conditions the 
synchrotron emission is characterized by a spectral index of the radio-to-optical 
continuum being approximately constant along the jet. A characteristic feature 
of the obtained broad-band synchrotron spectrum is an excess at X-ray frequen- 
cies, similar to the one observed in some objects by Chandra. As compared to 
the uniform jet models, the velocity shear across the radiating boundary region 
leads to decrease and frequency dependence of the observed jet-counterjet radio 
brightness asymmetry. We conclude that a careful investigation of the observa- 
tional data looking for the derived effects can allow to evaluate the role of the 
boundary layer acceleration processes and/or impose constraints for the physical 
parameters of such layers in large scale jets. 



Subject headings: acceleration of particles — galaxies: jets — radiation mechanisms: 
nonthermal 
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1. Introduction 

Jet radial stratification was proposed in order to interpret some radio and optical ob- 
servations of extragalactic jets. The existence of the jet boundary layer with a velocity shear 
was first suggested by the observations of jets in FR I sources. Owen et al. (1989) concluded 
that the jet morphology in M87 is dominated by the boundary layer more than by the shock 
structure. The constancy of the radio-to-optical spectral index along the jet and the similar- 
ity of radio and optical maps with emissivity peaking near the jet's surface strongly suggested 
the acceleration of radiating electrons in situ at the boundary region. Later polarimetry of 
M87 jet (Perlman et al. 1999) confirmed its complex structure, which contains knots with a 
transverse magnetic field dominating optical emission, and a shear layer which is bright in 
the radio band with a strong polarization suggesting a highly ordered magnetic field parallel 
to the jet's axis. Jet-counterjet surface brightness asymmetries and magnetic field structures 
of the other FR I radio galaxies, 3C 31 and 3C 296 (Laing 1996 and Hardcastle et al. 1997, 
respectively) , were interpreted in terms of a model in which the jets consisted of a relativistic 
spine surrounded by a cylindrical layer with a velocity shear. A study of optical counterparts 
of the radio jets in BL Lac objects PKS 2201+044, PKS 0521-365, and 3C371 (Scarpa et al. 
1999) emphasized the problem of the reacceleration of electrons radiating within the opti- 
cal band in the equipartition magnetic field, as the observed constant synchrotron spectral 
index along the jet could not be explained by the first-order Fermi shock acceleration. A 
map of the radio spectral index of another well known blazar, Mkn 501, revealed a boundary 
region of its jet with an unexpectedly fiat spectrum (Edwards et al. 2000). Once again, this 
suggested action of some reacceleration mechanism at the jet surface. 

Some FR II radio sources also reveal a spine - shear boundary layer jet morphology. 
For instance, radio observations of 3C353 (Swain et al. 1998) were interpreted in terms 
of a Doppler-hidden relativistic spine surrounded by the slower moving boundary layer. 
Polarimetry has shown no magnetic field component transverse to the jet axis in the boundary 
layer, in agreement with previous observations of radio linear polarization of quasars (e.g. 
Cawthorne et al. 1993). The magnetic field being parallel to the flow velocity is expected in 
the presence of a strong velocity shear (e.g., Kahn 1983). VLBI observations of another FR 
II object 1055+018 (Attridge et al. 1999) revealed a fragmentary but distinct boundary layer 
with the mean longitudinal magnetic field present where interaction with the surrounding 
matter takes place (i.e., where the jet bends). Optical and radio observations of the quasar 
3C 273 were interpreted in terms of a two component jet model, in which a fast-moving 
proper jet was surrounded by the slow-moving radio cocoon (Bahcall et al. 1995). Smooth 
changes of the radio-to-optical spectral index along the jet and a lack of correlation between 
the optical flux and the spectral index (Jester et al. 2001), strongly suggest reacceleration of 
the radiating electrons taking place within the whole jet body. Jester et al. also emphasized 
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the discrepancy between the light-travel time along the jet in 3C 273 and the lifetime of 
electrons responsible for its synchrotron optical emission. This discrepancy could not be 
removed by relativistic beaming or the sub-equipartition magnetic field. 

A difference between FR I and FR II radio sources is ascribed to respective jet power and 
the importance of its interaction with the surrounding medium (Ghisellini et al. 1998). In 
general, FR I jets are thought to decelerate and spread out on the tens-of-kpc scale with their 
apparent magnetic field dominated by the component made perpendicular due to propagating 
shocks (e.g., Laing 1996; Laing et al. 1999). On the contrary, jets in FR 11 sources seem to 
remain relativistic and well collimated even if they are far away from the nucleus. At the 
same time, their apparent parallel magnetic field configuration results from the interaction 
with the ambient medium, thereby creating the shear layer with an approximately constant 
(limited) thickness and dominating the observed polarization (Cawthorne et al. 1993; Swain 
et al. 1998). The relativistic velocity on a large scale leading to apparent superluminal 
motions is also present in intermediate objects hke M87 (Biretta et al. 1999). 

The broad band spectral properties of BL Lacs and FR I radio galaxies give us another 
argument in support of a jet stratification. By comparing the multiwavelength observations 
of these two classes of AGNs in a framework of a unification scheme, Chiaberge et al. (2000) 
found evidence for a significant boundary layer emission, dominating the radiative jet output 
in FR Is. It is because FR 1 nuclei are over-luminous as compared to the predictions of the 
unification scheme, and because the ratios of radio to optical luminosities in BL Lacs are 
inconsistent with the ones observed in FR I's. The inferred velocities of such boundary 
regions are significantly slower as compared to the velocities of the central spines, but still 
relativistic, in order to explain anizotropic emission observed from the FR Is' cores. 

One should note, that our knowledge of physical conditions at the relativistic jet bound- 
ary is still insufficient. Several types of viscosity for the relativistic jet with a shear layer 
which confines the beam due to viscous interaction with the ambient medium were inves- 
tigated by Baan (1980). Cosmic ray viscosity in a boundary layer with a velocity shear 
which influences the hydrodynamics of the jet and particle distribution along its surface 
was studied by Earl et al. (1988) and Jokipii et al. (1989). Analogous radiation viscosity 
was investigated by Arav and Begelman (1992). The issue of jet stability with respect to 
magnetohydrodynamic Kelvin- Helmholtz (KH) instabilities of a boundary layer in a vortex 
sheet approximation was studied by Ferrari et al. (1978, 1980, 1981); Hardee (1979, 1983) 
and Birkinshaw (1984), among others. The influence of the flnite thickness shear layer was 
considered by Ferrari et al. (1982) and Birkinshaw (1991). More recently, Hanasz and Sol 
(1996, 1998) investigated the growth of KH instabilities in two component relativistic and 
nonrelativistic jets. Such instabilities at the interface of a fluid beam and an external medium 
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can influence the apparent jet morphology (Ferrari et al. 1978; Hardee 1979; Benford et al. 

1980) and generate turbulent MHD waves which accelerate particles (Ferrari et al. 1979; 
Eilek 1982). Ferrari et al. (1979) demonstrated that a small amplitude Alfven turbulence 
developed by the KH instabilities in extended radio lobes can accelerate electrons to high 
energies via the wave-particle resonant interactions. Eilek (1982) showed that MHD tur- 
bulence generated by the surface instabilities and accelerating the electrons is restricted to 
a narrow turbulent edge of the jet. During the last few years, the numerical modeling of 
relativistic jets provided some additional information on the nature of such transition layers. 
3D simulations revealed the presence of a shear layer with a high specific internal energy and 
highly turbulent, subsonic and thin cocoon surrounding a relativistic flow (Aloy et al. 1999). 
Simulations of radio emission from such jets (Aloy et al. 2000), show a radial structure in 
both intensity and polarization. One should note that the reason for the generation of the 
turbulent shear layer in those simulations was numerical viscosity, which only qualitatively 
mimics the real viscous processes. 

Till the present time no one performed a more thorough study of both acceleration 
processes acting at jet boundaries and the resulting observational effects. Our present paper 
is intended to study the possibility that the turbulent boundary layers of large scale jets can 
substantially contribute to the radiative jet output. The modifications of the standard models 
of jet radiation arise from particle acceleration acting continuously within the considered 
region and kinematic effects due to the velocity shear. In the present study, we neglect a 
possible separate and different energetic electron population present within the jet spine. 
We concentrate on the boundary layer radiation, in order to evaluate a role of the electron 
acceleration processes in this layer and to constrain hardly known physical parameters of this 
region. Below, in Section 2 we discuss an acceleration mechanism creating at the jet boundary 
a power-law particle distribution with a much harder component at its high energy end. The 
radiation of such electrons through the synchrotron, synchrotron sclf-Compton and external 
Compton processes, modified by the kinematic eff'ects connected with the jet velocity shear, 
are derived in Section 3. Consequences of the considered acceleration mechanism in the 
context of multiwavelength large scale jets observations are discussed in Section 4. Finally, 
a few conclusions are presented in the last Section 5. 



2. Acceleration of cosmic ray electrons 

Particle acceleration in a velocity shear layer was proposed in the early eighties by 
Berezhko with collaborators and then gradually developed (c.f. a review by Berezhko 1990). 
The action of such a mechanism at relativistic shear layers occurring at side boundaries 
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of relativistic jets was discussed by Ostrowski (1990, 1998, 2000). The mechanism was 
considered to provide high energy cosmic rays and to be an important factor influencing 
the dynamics of relativistic jets in extragalactic radio sources, thereby forming cosmic rays 
cocoons. Ostrowski and Sikora (2001) proposed that such an acceleration mechanism could 
provide a cosmic ray proton population being a substantial pressure factor in FR II radio 
source lobes. Below we discuss the application of the model of Ostrowski (2000) to the 
electron acceleration in the large scale jets. 



2.1. Maximum electron energy 

Let us assume that a thickness of the transition layer between the considered tens- 
of-kiloparsecs scale jet and the surrounding medium, D, is comparable to the jet radius 
Rj ~ Ikpc (c.f. Owen et al. 1989; Swain et al. 1998, see also Aloy et al. 1999), and that 
the magnetic field B frozen into tenuous plasma at and outside the boundary, is parallel 
on average to the flow velocity, U. The scattering on the magnetic field irregularities in 
a turbulent medium with a perpendicular (= radial) mean velocity gradient results in the 
acceleration of energetic particles injected into the boundary layer. 

In our simple approach the efficiency of the acceleration process in the turbulent shear 
layer depends on the velocity structure of the layer and the particle mean free path A. In 
the considered case of the mean magnetic field aligned along the jet in a highly turbulent 
boundary layer, a particle gyroradius, r^, can be taken for A, and a relatively inefficient cross- 
field radial diffusion provides a particle escape mechanism from the acceleration region. For 
extremely high energy particles with Vg > D, the jet boundary can be approximated as a 
surface of a discontinuous velocity change. Such a particle moving near the jet boundary 
interacts with the magnetic field perturbations and ocassionally crosses the discontinuity 
thereby increasing its energy on average. Numerical simulations for 'typical' jet parameters 
(Ostrowski 1998) reveal that protons can reach energies above 10^^ eV in this mechanism. 
For particles with lower energies, < D, the acceleration within a finite thickness shear layer 
has to be considered. In this case, we deal with two different acceleration processes, the well- 
known second order Fermi acceleration in the turbulent medium and the acceleration due 
to cosmic ray viscosity in a velocity shear at the boundary. The cosmic ray viscosity in the 
context of particle acceleration was discussed by Berezhko in the early eighties (see Berezhko 
1990 and references therein) and then by Earl et al. (1988) and Jokipii et al. (1989) in the 
case of nonrelativistic flows, and by Ostrowski (2000) for relativistic flow velocities. Below, 
we consider these mechanisms to provide relativistic electrons radiating through synchrotron 
and inverse Compton processes. Because of rapid energy radiation losses at high energies. 
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the electrons are always expected to satisfy the required condition Vg <^ D. Then, the 
acceleration time scale in a highly turbulent shear layer can be estimated (Ostrowski 2000) 
as 

T,^ ^ ( !a) (2.1) 

where y is a characteristic velocity of the magnetic turbulence, U is the jet velocity and we 
put Tg for A. The first term in the denominator represents the ordinary Fermi mechanism 
and the second one represents the viscous acceleration. The Fermi process dominates when 
the particle energy is small enough to satisfy Vg < D{V/U). 

We consider the jet velocity to be comparable to that of the light velocity, c, even on large 
scales. The velocity V is comparable to the Alfven velocity, Va, for a subsonic turbulence and 
possibly a few orders of magnitude smaller than U. MHD turbulent modes in a shear layer 
can be dominated by Alfven waves, because the other types of waves are damped much more 
rapidly. In order to estimate the Alfven velocity, Va — B / y/^np, we assume that the mass 
density p of the jet and its extended lobes is dominated by nonrelativistic protons, p = iripn 
(Sikora and Madejski 2000). With the anticipated parameters for the boundary layer of the 
tens-of-kpc scale jet - -B ~ 10^^ G and the proton number density n ^ lO^'^cm"'^ (c.f. 
Ferrari et al. 1979) - one gets Va ~ 2 ■ lO^cm/s. Then, with the condition Vg < D(y/U) 
satisfied and V = Va, the characteristic acceleration time scale (2.1) for electrons ('e') reads 
as 2 

Tacce - ^-^ {^) ^5-10'jB;^Vi' [s], (2.2) 

where 7 is the electron Lorentz factor, Vg = V^/lO^cms^^ and Bf^c is the magnetic field 
in micro-Gauss units. Instead of n, we will use Va to parametrize the acceleration process, 
because the mass density is not an observed quantity. One should note, that the equation 
(2.2) introduces an optimistic scenario for the acceleration process involving the large am- 
phtude MHD turbulence at the scale of r^, leading to short scattering free path A ~ r^. The 
time scale for particle escape from the acceleration region due to cross-field diffusion can be 
evaluated as : ^ 

7;.c ~ — = 6 • 10^3 7"' B,G [s], (2.3) 

K,± 

where k± is the cross-field diffusion coefficient, k± = rjrgc/S, with rj < 1 being a numerical 
scaling factor, and D is the acceleration region (boundary layer) thickness {D^pc = D/1 kpc). 
Particle escape determines the electron energy spectrum, when the electrons satisfy the 
condition Tacc,e ~ T'esc , i-e., when Vg {3/riy^'^D(V/c). For the considered boundary layer 
parameters and a strong turbulence condition 1, radiative cooling becomes important 
at much lower electron energies. We consider radiative losses due to the synchrotron and 
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the inverse-Compton radiation (in the Thomson regime), to yield the loss time scale 

where X is a ratio of the photon field energy density to the magnetic field energy density. 
Therefore, the maximum electron energy in the acceleration process can be estimated by 
comparing the time scale (2.2) with the scale for radiative losses, (2.4). For a value of 
X <1, possibly vahd at the considered large distance from the active galactic nucleus (see 
next section) , the maximum electron energy given by the condition Tacc,e — Tioss is 

7e, ~ 4 • 10« B;^/^ (2.5) 

In our model, a rapid radiative cooling of electrons with such high maximum energy is 
compensated by acceleration taking place along the jet within the whole boundary layer 
volume. For lower turbulence amplitude, the acceleration process can proceed slower to yield 
a lower value for 7eg. One may note, that at 7 = ^^.q the ratio Tf.sc/Tacc ~ 7 • IQ^D^^p^B^Q 
(~ 7 • 10^ for the considered jet parameters). 



2.2. Electron energy distribution 

A simple time dependent model of Ostrowski (2000) uses a regular acceleration term 
to represent continues second-order Fermi acceleration plus the not-considered here oblique 
shocks possibly formed at the jet side boundary. After switching on the injection of seed 
electrons with the low initial Lorentz factor 70, a spectral evolution of the distribution 
function ne{j,t) for particles accelerated in the shear layer begins as a power-law with a 
growing energy cut-off. Then, either it forms a stable high energy cut-off when particle escape 
becomes substantial at some 7 < jeq, or it evolves into a characteristic shape consisting of two 
different components: a power- law part at low energies, rig (7) oc 7"'^, finished with a hard 
component modeled here as a nearly monoenergetic peak at 7 = 7eg due to the accelerated 
particles' pile-up caused by losses (Fig. 1). Normalization of the electron energy distribution 
depends on an unknown electron injection rate and is a free parameter of the model. In 
principle, the particle high energy component can grow to large values, limited only by the 
inefficient side way particle diffusion. However, when the energetic particle pressure becomes 
comparable to the ambient medium pressure, it can modify the transition layer structure. 
Then, some non-diffusive cosmic ray transport (turbulent diffusion, convective motions) or 
other non-linearities of the process can stabilize the average spectral distribution at the 
boundary on long time scales, or result in its locally intermittent character. Below, because 
of the very nature of the acceleration process, we consider the stationary avarage electron 
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spectrum formed within the boundary layer consisting of two basic components: a flat power- 
law component at lower energies and a high energy bump. Using the Haeviside function, 0, 
and the Dirac delta function, S, at energies above the injection energy, 7 > 70, the considered 
electron distribution (cf. Fig. 1) can be expressed by an analytic formula 

ne(7) = a7-" Qile, - 7) + &5(7 - 7e,), (2.6) 

where a and b are normalization parameters, and we put a = 2. In the following discussion, 
we use numerical values for a and b providing equipartition between the magnetic field energy 
density ub = B'^/Stt and the energy densities of both electron spectral components: 

/ (7mc^)a7~^e(7e5-7)d7= / {-fmc^)bS{-f - -feq) d-f ^ -ub, (2.7) 

J JO J 10 

where ^eq > 7o- With assumed B — 10~^ G and maximum electron energy much higher than 
the injection energy, •jeq ~ 10^ 3> 70 ~ 10 — 100, the equipartitional normalization of the 
power-law electron distribution and the number density of the monoenergetic peak electrons 
are a ~ 10~^cm~^ and b ~ 10~^"'cm~^, respectively. 

Let us mention, that an analogous model was discussed by Schlickeiser (1984), who con- 
sidered similar pile-up mechanism for the relativistic electron distribution and studied evo- 
lution of the formed monoenergetic electron peak by solving stationary and time-dependent 
diffusion equation in the momentum space. In his model, the first-order shock acceleration 
accompanied by the second-order turbulent Fermi acceleration creates the pile-up growing 
bump at the energy, where the radiative losses balance the energy gains. Schlickeiser assumed 
the mean free path of the electron being independent on its energy, contrary to the consid- 
ered by us A oc 7. It was shown, that the stationary two-component electron spectrum (the 
power-law plus the monoenergetic peak) forms in a presence of efficient shock acceleration, 
when the escape time scale is much longer than the Fermi acceleration scale. After switching 
off the shock acceleration, at the late phase of the evolution the high energy electron peak 
evolves into the ultrarelativistic Maxwell-Boltzmann distribution with the absolute height 
of the bump decreasing as exp(— t/Te^c)- 

Below, we assume that the (quasi-) stationary averaged spectrum of electrons formed at 
the jet boundary layer can be represented by the distribution (2.6). As long as the diffusive 
particle escape to the sides is inefficient, the discussed acceleration process will create a high 
energy particle pile-up in a natural way (cf. Stawarz and Ostrowski 2001). The distributions 
produced within the jet boundary layer can vary in time and space reflecting the changes of 
the local conditions due to the jet interaction with the surrounding medium and action of 
non-linear effects in the acceleration process. However, the form of the averaged observed 
electron spectrum must consist of the two aforementioned components: the more or less exact 
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power-law and the broadened by the fluctuating background high energy bump preceding 
the distribution cut-off. Without a detailed acceleration model available, depending also on 
the unknown velocity structure of the boundary layer, we use the simple form (2.6) to discuss 
the role of these two spectral components for a particular example of comparable energies 
stored in both components. 



3. Radiation of shear layer electrons 

Wc consider two radiation processes taking place at the relativistic jet boundary - the 
synchrotron (SYN) emission of ultrarelativistic electrons with the energy spectrum (2.6), 
and their inverse Compton (IC) cooling. In the case of magnetic bremsstrahlung, we assume 
randomly orientated magnetic fields within the turbulent medium of the transition layer. 
A synchrotron self-absorption is neglected (i.e., the assumed optical depth for this process 
Tgyn = Hsynl ^ 1; whcrc Hsyn IS a syuchrotrotrou self- absorption coefficient and I is an 
emitting region size). Therefore, in our calculations we limit the considered frequencies 
somewhat arbitrarily to the 'safe' range v > Vabs — 10^*^ Hz, where self- absorption effects 
are unhkely to occur in the considered large scale jets. Additionally, the present analytic IC 
scattering calculations are limited to the Thomson regime, assuming a step scattering cross 
section a = for 7 ej < 1 and otherwise, where is the seed photon energy expressed in 
units of the electron rest energy, m(? . A high energy spectral cut-off due to absorption of 
gamma rays via photon-photon pair creation on Cosmic Infrared Background (CIB) photons 
is expected to occur for distant sources (e.g., Renault et al. 2001). As we do not consider 
such processes, the obtained spectral distributions above TeV energies are the ones derived 
for close vicinity of the radiating jets and cannot be directly compared to observations. 
Considering the source of the target photons in a tens-of-kpc (or larger) scale jets, we neglect 
a possibility of comptonization of the galactic narrow line radiation and IR dust emission 
(Celotti et al. 2001), and we study the synchrotron self-Compton (SSC) radiation and the 
Compton scattering of external CMB photons (EC). With the assumed B = 10~^ G, the 
magnetic field energy density at the jet boundary is Ub = -B^/Svr ^ 4 ■ 10~^^ erg/cm^. The 
energy density of the CMB radiation in a source frame moving with a bulk Lorentz factor 
r (neglecting the cosmological redshift correction) is Ucmb = ^cmft^^ ~ 4 • 10"-*^^ erg/cm^, 
where an asterisk denotes quantities in the observer's rest frame. In the jet rest frame, 
Ucmb can be of the same order of magnitude as ub, and both processes can be significant 
even for moderate F. One should note that in such a case, the estimate (2.5) for still 
holds. Celotti et al. (2001) considered also comptonization of the blazar ('bl') emission by 
electrons in the boundary region of the kpc-scale conical jet, with a low average boundary 
layer Lorentz factor F. The energy density of such radiation as seen in the boundary layer 
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local rest frame is uu = F^^ / ^ttc fa 3 • 10~^° z^^^ F"^ erg/cm^, where L^^ ~ lO^^erg/s 
is the intrinsic synchrotron blazar luminosity, Tp^ ~ 10 is the bulk Lorentz factor of the 
parsec-scale jet, and 2; is a distance from the nucleus {z^pc = z/kpc). The blazar emission 
dominates over the magnetic field energy density at distances z^p^ < lO/F, and over the 
CMB field at the scales Zkpc < SO/F^. Hence, for the studied large scale jets {z > lOkpc) 
and the relativistic fiows (F > 2) the blazar emission is neglected. In order to derive analytic 
formulae for emission we are forced to make a simplifying approximation of the isotropic 
distribution of the synchrotron radiation in the source local rest frame. In reality, the 
shear layer synchrotron radiation distribution is anizotropic due to kinematic effects of the 
nonuniform relativistic flow. A numerical modeling of such anizotropic radiation field will 
be presented in a forthcoming paper. 



3.1. Synchrotron radiation 

The synchrotron emissivity of an isotropic ultrarelativistic electron distribution Ueij) 
averaged over a randomly orientated magnetic field, B, is 



Jsyn[l^) = — / R 2 —A (3-1 



where ci = 3eB / Anmc and R{x) is given by a combination of the modified second order 
Bessel functions (Crusius and Schlickeiser 1986) as 



^2 



fl(^) - \ if./, (I) K,,^ (I) - 0.3 \ [Kl,, (I) - A1,3 (|)] • (3.2) 

Upon integrating the equation (3.1) for the electron distribution (2.6) and expanding the 
function R(xeq) for Xf.q = T^/cillq ^ 1' R{^eq) ~ 1-8 " ^Iq^ , One Can obtain simple analytic 
expression for the synchrotron emissivity. For x^q ^ 1, the function (3.2) decreases expo- 
nentially, R{xeq) oc exp(— Xeg), and thus we neglect emission at frequencies u S> cij"^^. For 
parameters a_7 = a/(10~^cm~^) and 6_i4 = 6/(10~^^cm~^), the synchrotron emissivity is 

Jsyniiy) = [fl y-^''' + /2 e(z/,,„,eg " v). (3.3) 

where 

/i = 3 • 10-36 Bjg a_7 [cgs], 
/. = 2.10-^^-^/'<j6_,4 [cgs], 

Vsyn,eq^^-'^B^,Gllq [Hz] . 
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For illustration, let us compare the synchrotron radiation spectra of the particle distributions 
in the acceleration model of Ostrowski (2000) in successive times from the switching on of 
the acceleration process (Fig. 1). From the beginning of injection, the accelerated electrons 
form a power-law spectrum with an energy cut-off which grows with time. Then the resulting 
emission has (for a = 2) a power-law form oc z/^^/^ with a growing cut-off frequency. When 
the maximum electron energy approaches 7eg, a particle density peak starts to form and 
grow at 7 ~ 7eg. Then, an additional radiation component oc v^l^ appears and eventually 
starts to dominate the spectrum at highest synchrotron frequencies, reaching the maximum 
at Vsyn,eq (~ 10^^ Hz for 10^ and ~ 10). 



3.2. Synchrotron Self-Compton radiation 

The inverse- Compton photon emissivity in the source frame, hic{e^ fi), is given by 

na., Si) = c / *. / ia„ I d, jdilA^-p COS i,) a n*., Si.) il^), (3.4) 

where e is the photon energy in the electron rest mass units, e = hiy/mc^, ni{ei, Qj) is the 
seed photon number density (wc put the indices i = syn, ic — ssc for SSC emission, and 
i — cmb, ic — ec for EC emission), ne(ee, l^g) is the electron energy distribution, and ip is 
an angle between the electron and the incident photon directions (e.g., Dermer 1995). All 
quantities are given in the source rest frame, which in our case, is a part of the boundary 
region - a cylindrical layer moving with a constant Lorentz factor F (see section 3.4). The 
rest-frame emissivity in cgs units can be found next as jidy, ^) = hehic{e, Q). 

Computations of the SSC emission arc performed in the Thomson regime with a scat- 
tering cross section independent of the seed photon energy, and with the scattered pho- 
tons beamed along the electron direction. The synchrotron radiation and the electron 
distribution are assumed to be isotropic in the source frame, ne(7,f2e) = ^e(7)/47r and 
''T'syni^sym^syn) — f^syni'^syn) / ■ Thus, the average energy of the comptonized synchrotron 
photon is e = ^'y'^^syn, and in equation (3.4) one can put a — ut ^iS^ — (e — \'i^^syn) to 
obtain 

nssc{^,^ = ^ y desyn j dj Usyniesyn) rieij) S (^6- ^-/"^esyr^ ■ (3.5) 

For the electron energy distribution given by the equation (2.6) and the optically thin syn- 
chrotron intensity, Isyn{i^syn) = jsyn{j^syn) I = ^ syn n syn{(^ syn) (whcrc / is an emitting region 
linear size), putting In (z/^aa;/i^a&s) ~ 10, where v^ax = niin(i/s2^„,eq , mc^/h'yeq) = mc^/h'yeq, 
and neglecting absorption due to photon-photon pair production, one can obtain the isotropic 
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SSC emissivity in the Thomson regime 

Jsscil^) = [/3 1^-'^' + /4 1^'^'] e{l^ssc,eq ' 1^) , (3.6) 

where 

/3 = 10-^^ B^J^ (3 + 10-^ 7e, a_7 6_i4) [cgs] , 

/4 = 5 • 10-^« B'J^ 7;//-^ [cgs] , 

4 

i^5sc,eq = g J^max 1% = 1.6-10^° 7eg [Hz] . 

Note, that the high energy electron bump contributes to both components of the SSC emis- 
sion (i.e., also to the one oc v"^^"^). 



3.3. External Compton radiation 

Cosmic microwave background (CMB) intensity, Icmb{^), has a blackbody form with a 
temperature T* = 2.73K. In a source frame moving with the relativistic Lorentz factor F, 
the CMB radiation is strongly anizotropic, providing an external radiation flux mainly from 
the direction opposite to the jet direction. Hence, for ultrarelativistic electrons, one can put 
(1 — /? cosV') = (l+z^e), where fig = {cos~^ 0e)- Then, the scattered photon energy is e = 
7^ ^cmb (1 + fJ'e)- Therefore, in the Thomson regime a = ax — Vt^.) 5[e — ^"^e^mb (1 + A'e)]) 
and for the isotropic electron distribution, the equation (3.4) reads as 

fiecie, n) = ^{1 + ii) J decmb J dj ricmbi^cmb) ne{j) S[e- 7^ec,n6(l + n)] . (3.7) 

Considering the power-law component of the electron energy distribution, it is convenient to 
model the CMB field as a monochromatic radiation with photon energies €cmb = T < >= 
2.7 • r A; r*/mc^, and a number density ncm6 = rti^^/(< > mc^). For the monoenergetic 
electron peak component, one should use the exact blackbody spectrum transformed to the 
source frame. Hence, the emissivity of the comptonized CMB radiation, jec{j^, is 

€)(l^ec,eq - i^), (3.8) 

where 

/5(/i) = 8-10-^ia_7[r(l + ^)]^^/2 [cgs], 

/6(/x) = io-«S;/^-i4(i + /^)-' [cgs], 
M/.) = 3-io-^^r(i + /.)-S;/ [Hz-^j, 



exp[/7(/^) H - 1. 
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iyecM^2-10''r{l + f,)^l^ [Hz]. 

The high energy bump radiation with frequencies u < fT~^{iJ>) (i.e., anizotropic Rayleigh- 
Jeans part of the blackbody CMB spectrum comptonized by the monoenergetic electrons) 
has a power-law form oc i/^. 

3.4. Kinematic effects 

Let us briefly discuss the main imphcations of the kinematic effects on the jet radiative 
output. First of all, the shear layer consists of sub-layers moving with different velocities and 
thus, generating different beaming patterns. In order to specify the jet velocity structure, 
we assume a uniform flow Lorentz factor Tj inside the jet spine, i.e. for < r < Rj where r 
is a distance from the jet axis, and a radial profile F = F(r) within the cylindrical transition 
region Rj < r < R^, where the flow velocity decreases and becomes zero (= the ambient 
medium velocity) at the radius Rc- Let us consider the simplest case of the linear dependence 

F(r) = l + (F,-1)AZ^. (3.9) 
it(, — itj 

for Rj < r < Rc- In each considered sub-layer, in the 'local source frame' the synchrotron, 
SSC and EC emissivities are approximately given by formulas (3.3), (3.6), and (3.8), respec- 
tively. Assuming that the radiating electrons are distributed uniformly (as seen in the local 
rest frames) inside the whole boundary region, one can find the observed emissivity of the 
sub- layer at the distance r from the jet axis as 

,.•(... r.r, ^ .-(nr),. . - . (3.io) 

where the index i — syn, ssc or ec, and S{r, 6*) is a Doppler factor for a sub- layer moving 
with the Lorentz factor F(r) = [1 — /9^(r)]~^^^ at an angle 9* = cos~^ //* with respect to 
the line of sight, S{r,6*) = 1 /F(r) [1 — j3{r) fi*]. Next, the observed flux density from the 
boundary layer volume V* can be found as S*{i'*,6*) = J dV* j*{i'*,9*,r), where d is 
the distance to the observer. Noting that F(r)[l + fi] = d[l + /U*]/[l + l3{r)], and neglecting 
the slowly varying factor [1 + /i*]/[l + (3{r)] (c.f. Dermcr 1995) from the above equations one 
obtains the flow modified beaming patterns for the synchrotron, SSC, and EC radiation, as 

S:y^{u*,9*) cx (i.*)-° r drr5^+-{rX), (3.11a) 

rRc 

S*,,{u*,9*) <x{iy*y drr S^+"{r, 9*) / sin 9*, (3.11b) 
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fRc 

Sl^{p\e*) (X drr6'^+^'^{r,e*), (3.11c) 

jRj 

where a is the spectral index of the respective power-law emission, and a — —2 for the 
Rayleigh- Jeans part of the comptonized CMB spectrum. In the case of SSC emissivity, we 
assumed isotropic distribution of the seed photons in the source rest frame, and we put the 
effective emitting region linear size scaling roughly like I cc D / sin^^*. Integration in (3.11) 
can be performed numerically for the given values of Tj, Rj, R^, 6*, and for the assumed 
radial velocity profile. 



3.5. Composite spectra 

Our model of high energy radiation generated within the jet boundary region has a few 
free parameters. They include normalization of the power-law electron distribution, a, a 
height of the monoenergetic electron peak, b, the magnetic field induction, B, the Alfven 
velocity in the shear layer, Va (or the mass density of the plasma within the boundary region, 
p), and the injection electron energy, 70. In general, the electron power- law spectral index, 
(7, can also be treated as a free parameter represented in this paper by the model value 
a — 2. Here, Va (or p) and B limit a value of the maximum electron energy, ^eq- Relative 
normalizations of radiation spectral components - synchrotron, SSC and EC emission - 
depend on the above parameters and hence the observed composite boundary layer spectral 
energy distribution varies with the selected values of a, b, B, j^q, 7o together with the 
geometrical factors Tj, Rj, Rc and 9*. 

For a 'typical' tens-of-kpc scale jet the magnetic field and the Alfven speed are, say, 
B ^ 10^^ G and Va ~ 10*^cm/s, and the maximum electron energy can be as high as 
7eg ~ 10*^ (Eq. 2.5). We assume, that in the boundary region 70 ^ 7eg and that the 
electrons are uniformly distributed within the whole radiating volume (as seen in the local 
layer rest frame) . The magnetic field energy density can provide estimates for normalizations 
of both electron spectral components, if we require the equipartition: a ~ 10~^cm~^ and 
b 10~^^cm~^ (Eq. 2.7). The thickness of the large scale jet boundary layer and the jet 
radius are assumed to be of the order of a kiloparsec, D ^ Rj ^1 kpc, and hence Rc — 2 Rj. 
Below, for illustration, we consider a part of the jet with the observed length « 1 kpc. 
We also assume, that the jet spine is highly relativistic even at large scales, with that the 
bulk Lorentz factor Tj ^ 10 (c.f. Ghisellini and Celotti 2001). The jet with the described 
parameters is used to illustrate possible role of kinematic effects at the observed brightness 
distribution along the jet. 

The observed spectral energy distribution (SED) of the emission generated at the con- 
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sidered part of the boundary region, i/*S*{i'*,9*), is shown on figure 2 for four different 
viewing angles. We put the source distance from the observer d ~ lO^kpc. In order to 
estimate the observed luminosity of the boundary layer emission, one has to integrate over 
the intrinsic luminosities of the sub-layers with volumes dV* = 27i Az r dr and the Doppler 
factors d{r, 9*). For example, the synchrotron luminosity of the boundary layer electrons is 



where 7ioss = l/Tioss denotes the synchrotron energy losses. In our illustrative example, 
synchrotron luminosity of the high energy bump electrons observed at 9* — 45° (i.e. the 
observed flux multiphed by the Ai^d"^) can be as high as lO^^erg/s. 



In this section we compare the model of a radiative boundary layer with the large scale 
jet observations. We do not consider the possible contribution from the separate population 
of relativistic electrons accelerated within the knots due to the first order shock acceleration. 
Our considerations are intended to describe only the main spectral characteristics of the 
relativistic electrons accelerated within the turbulent boundary region. 



For the anticipated large scale jet parameters, assuming equipartition between the mag- 
netic field and the relativistic electrons accelerated within the shear layer, one will find the 
power-law electron component dominating the radio-to-optical jet boundary emission, as il- 
lustrated on figure 2. As a result of a continuous acceleration taking place within the whole 
boundary layer, the radio-to-optical spectral index is constant along the jet. Of course, 
changes of the physical conditions with the distance from the nucleus (e.g., changes of the 
magnetic field, turbulent conditions, etc.) can result in smooth changes of the slope of the 
synchrotron continuum, hke in the case of M 87 or 3C 273 (Meisenheimer et al. 1996; Jester 
et al. 2001, respectively). In our model the acceleration process creates the high energy 
electron spectral bump and the optical emission joins smoothly the X-ray bump emission. 
Such spectral character is different from the case usually considered in the literature of the 
power-law continuum ended by the cut-off at optical or X-ray frequencies, followed by the 
inverse- Compton component. 




syn — 




(3.12) 



4. Multiwavelength large scale jet boundary emission 



4.1. Radio-to-optical continuum 



The important manifestation of the boundary layer radio-to-optical emission is a de- 
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crease of the jet-counter jet brightness asymmetry as compared to the uniform jet, suggested 
previously by Komissarov (1990). This effect arise from the fact, that in a presence of the 
velocity shear, for a given viewing angle 9* > l/Tj the observed synchrotron emissivity is 
maximized for r(r) = 1/ sin^* (i.e., for the maximum value of the Doppler factor 6{r,6*)), 
while the uniform jet radiation, if present, is Doppler-hidden. Figure 3 illustrates the bound- 
ary layer jet-counterjet brightness asymmetry, 

S* (0*) J^'drrS^+^{r,0*) 

for the spectral index a = 0.5 corresponding to the radio emission of our power-law electron 
component and the Lorentz factor profile (3.9) with = 2Rj and Tj = 10 or 3. For 
comparison, respective fiux asymmetries for the uniform jet models are also plotted. Due to 
the considered effect, for = 10 and the moderate inchnations, the observed jet asymmetry 
can be reduced by more than one order of magnitude in comparison to the uniform jets. 
Hence, the middly-relativistic velocities inferred from the observed jet asymmetries of the 
tens-of-kpc scale jets can correspond to the slower boundary layer and not necessarily to the 
fast spine, which can still be highly relativistic at the observed distances. 



4.2. X-ray bump emission 

In the model presented here, the X-ray radiation of large scale jets is dominated by the 
high energy electron bump spectral component. As illustrated on figure 2, the observed 
synchrotron X-ray fiux can be above the extrapolated radio-to-optical continuum. The 
spectrum inclination at X-ray frequencies can be different from the power-law observed at 
lower frequencies, and the difference depends also on a viewing angle. Figure 4 illustrates this 
dependence for the effective X-ray spectral index ax,eff{d*) computed between hu^ = IkeV 
and /t = 5 keV as 

OI.X,eff[9 ) = \ — r ^ / ^1 ■ (4.2) 

where S*x{i^*,0*) oc / r dr S'^{r,6*) R{xlg), a function R{xl^ is given by equation 3.2, and 
x*q — V* /ci'y1q5{r,9*). The difference between oix,eff for different jet Lorentz factors and 
different viewing angles results from the spectral character of the high energy electron bump 
emission, which is strongly concentrated around the maximum synchrotron frequency Vgyceq- 
Transformation to the observer's rest frame, which depends on the flow Lorentz factor and on 
the inclination 6*, results in the fact, that for the fixed observed frequency u* in this spectral 
range {v* ~ i'syn,eq ) different parts of the X-ray continuum can contribute to the radiative 



-17- 



output. For the higher Doppler beaming factors (i.e., higher Fj or smaller 9*) the effective 
X-ray spectral index of the bump electrons with ^ygg ~ 10^ approach the asymptotic value 
—1/3. In the particular case considered by us the counterjet has a steeper X-ray continuum, 
as compared to the jet spectrum. Moreover, for large viewing angles. X-ray emission of the 
jet softens significantly, opposite to the counterjet emission. Note, that such a behavior is 
opposite to what we expect in the case of IC models involving the low-energy tail of the 
electron distribution, in which the X-ray spectral index of the jet and the counterjet should 
be the same, and approximately constant with the distance from the nucleus and with the 
viewing angle. The boundary X-ray emission at the large scales can therefore vary among 
different types of the jets, reflecting different physical conditions at their edges and different 
kinematic effects involved. Thus, it is interesting to note that the large scale X-ray emission 
observed from several radio loud AGNs exhibit the different spectral characteristics. 

Chandra X-ray Observatory has detected significant emission from the large scale jets 
in quasars (e.g. PKS 0637-752: Chartas et al. 2000), FR II sources (e.g. Pictor A: Wilson 
et al. 2001), and FR I radio galaxies (e.g. M 87: Marshall et al. 2001b). In the majority 
of cases, thermal bramstrahlung can be excluded, and the non-thermal models seem to 
provide the only possible explanation for the X-ray emission (Harris 2001). However, the 
spectral index of the X-ray non-thermal continuum is different for different sources. Also 
the ratio of the X-ray flux to the optical flux changes signiflcantly in the Chandra jets 
sample. This suggest, that different radiative processes can be at work. The models of 
the non-thermal large scale jet X-ray emission usually involve the synchrotron emission of 
high energy electrons, with Lorentz factor 7 ~ 10^ or the inverse Compton scattering of 
CMB photons by much less energetic electrons (7 ~ 10^) in a rclativistic jet. Both classes 
of models are subject to difficulties in explaining the variety of Chandra observations. The 
synchrotron emission from one population of the electrons accelerated in a rclativistic shock 
is excluded if the X-ray flux is above the extrapolated radio-to-optical continuum. Also, 
the continuous along the jet character of the large scale X-ray emission is inconsistent with 
short lifetimes of the electrons radiating X-rays in the equipartition magnetic flelds, unless we 
consider continuous acceleration within the whole jet body. On the other hand, the 'beamed' 
external Compton (EC (&eam)) models require highly rclativistic jet flows on the tens-of-kpc 
scales, and extrapolation of the single power-law electron energy distribution to low energies, 
while, in most of the cases, there is no direct evidence supporting these assumptions. As 
discussed in detail by Harris and Krawczynski (2001) such models can be excluded if the 
X-ray spectrum is much steeper than the one observed at radio frequencies. One should also 
note, that the long life-time of the electrons scattering CMB photons is in contrast with the 
knotty morphology of the X-ray jets, as well as with the off-sets between the radio and the 
X-ray peaks reported in several brightest knots. 
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The EC{beam) model with significant Doppler beaming seems to work well in the case 
of the quasar PKS 0637-752 (Tavecchio et al. 2000), where VLB! observations suggest the 
parsec-scale flow Lorentz factor Tj > 18 and the inclination angle 9* < 6.4°. For other 
quasars the situation is more complicated. For instance, observations of the jet in 3C 273 
(Marshall et al. 2001a; Sambruna et al. 2001) are insufficient to choose between the syn- 
chrotron or the EC origin of the X-ray knots' emission, and the existence of an additional 
electron population radiating in this range is not excluded. One should also note recently 
pubhshed HST observations (Jester et al. 2002) showing smooth transition from optical 
through UV to X-ray frequencies, as expected in our model. The X-ray luminosities of the 
known quasars' jets are in the range 10^^ — 10^^ erg/s, and their X-ray spectra are very flat, 
i.e. ax ~ 0.23 for 3C 207, (Brunetti et al. 2001), ~ 0.5 for PKS 1127, (Siemiginowska et al. 
2002), or ~ 0.8 for 3C 273 and PKS 0637. In our model, the observed luminosity of the X-ray 
emission generated at the boundary layer can reach the observed values for reasonable jet 
parameters (c.f. discussion in section 3.5), and the low effective spectral index is consistent 
with small jet inclinations (c.f. Fig. 4). The significant X-ray emission is also observed from 
the large scale jets in radio galaxies. The lower X-ray luminosities of such jets, lying in the 
range 10^^ — 10^^ erg/s, and their relatively steep X-ray spectra, ax ~ 1-0 — 1.5 (Kraft et 
al. 2000 for Cen A, Hardcastle et al. 2001 for 3C 66B, Worrall et al. 2001 for B2 0206 and 
B2 0755), correspond to the large inclination angles on figure 4 and/or lower jet Lorentz 
factors. In most of jets in the radio galaxies, the EC{beam) models are excluded, as they 
require too high beaming factors, the magnetic fields much below the equipartition values 
and too steep low-energy electron spectra. On the other hand, the synchrotron radiation of 
the shock accelerated electrons cannot explain the strong diffusive emission extending few 
kpc from the nucleus in, e.g., Cen A and 3C 66B. 

We expect, that the detailed investigation of the knots emission can allow verification 
of our two-component electron spectrum model, if the knots form by the shock compression 
of the energetic electrons previously accelerated at the boundary layer upstream of the shock 
(c.f. Begelman and Kirk 1990). More directly, the boundary layer X-ray radiation should 
manifests at the inter-knot sections of the jet. Due to the acceleration taking place con- 
tinuously within the whole radiating volume, such an X-ray flux should be approximately 
constant (or change slowly) along the jet. Unfortunately, there are no observations in large 
scale jets of the counterjet X-ray spectra. In some cases, the lower limits for the X-ray 
brightness asymmetry can be estimated, like for the 3C 66B (> 25, Hardcastle et al. 2001), 
Pictor A (> 15, Wilson et al. 2001), or PKS 1127 (> 5, Siemiginowska et al. 2002). There 
is evidence of a non-zero X-ray flux from the counterjet in the last object. Note, that in a 
framework of radiative boundary layer model, the jet-counterjet X-ray brightness asymmetry 
is smaller than the one in EC models with strong beaming. 
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4.3. Very high energy 7-ray emission 

Very high energy (VHE) 7-rays are directly observed only from a few nearby AGNs 
during their flaring state (c.f. a review by Kifune 2001). Very short time scales of the 
spectral variations at TeV photon energies suggest small sub-parsec sizes of the emitting 
regions, where the violent particle shock acceleration take place. In our model, the electrons 
accelerated within the boundary layer of the large scale jets are also able to produce such 
energetic radiation, mainly due to the IC scattering of the CMB photons, as illustrated 
on figure 2. One shoTiId note, that figure 2 corresponds to the most optimistic scenario, 
with the highly relativistic jet spine and very efficient acceleration creating the electrons 
with large Lorentz factors up to ~ 10^. Smaller values of Tj and would result in 
decreasing the observed VHE flux and shifting the spectral cut-off to the lower frequencies. 
Alternatively, for the lower jet velocities or smaller jeq at the few-of-kpc scales the blazar 
emission illuminating the jet from behind can provide seed photons to produce TeV radiation. 
Yet another possibihty considered in the hterature is producing the VHE 7-rays by the 
relativistic protons accelerated in the large scale relativistic jets (Aharonian 2002). 

VHE 7-rays generated at the distant jets lose their energies due to absorption on the 
GIB radiation, GMB photons, and other photon fields surrounding the source. Aharonian 
et al. (1994) suggested, that the electromagnetic cascades initiated by such an absorption 
create the extended pair halo around the 10-100 TeV extragalactic sources. VHE radiation 
generated at the jet boundary layer should naturally contribute to this process. 

4.4. Discussion of the free parameters of the model 

In order to illustrate the main physical features of the considered model we selected 
its parameters in a way to exhibit clearly its main characteristics. In particular we neglect 
possible radiation components from the jet spine and the ambient medium next to the 
boundary layer, as well as from the shocks possibly formed in the fiow. In real jets the 
respective components add to the observed object radiation. For the considered radiating 
electrons we assume energy equipartition with the background magnetic field and the equal 
division of energy between the power-law and the final bump components. One should note 
that application of this simple two-component distribution allows for analytic study of the 
formed spectrum, it enables a clear separation of signatures of these both components in 
it, but to model the observations a more elaborated distributions can be required. Such 
distributions will involve an initial power-law segment in the spectrum, which flattens at 
larger energies and forms a flnal (not sharp) bump before the exponential cut-off (work in 
preparation; cf. Petrosian and Donaghy 1999). We also take a particular set of background 
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conditions for the jet and its boundary layer (the magnetic field and the matter content, the 
turbulence and the large scale kinematic structure) and the ambient radiation fields for this 
model. Let us consider how the main predictions of the presented model depend on certain 
values of the free parameters assumed in the text, and in particular if the maximum energy 
of the boundary layer electrons can substantially differ among different jet models. 

For the efficient particle acceleration within the shear boundary layer we assume its tur- 
bulent character with high amplitudes MHD or plasma waves scattering energetic electrons. 
The assumed configuration of the mean magnetic field within the shear region with dominat- 
ing components parallel to the jet axis was already discussed and seems to be confirmed by 
the observations. For these conditions the acceleration process acts as described in Section 
2, with the time scales for particle escape and radiative losses (Eq-s 2.3 and 2.4) satisfying 
the inequality T^gc > Tioss, equivalent to D^pc > O.OAB'q^. One may note that for a range 
of magnetic fields in large scale jets, 10 < B^q < 10^, the above condition can be easily 
satisfied for the transition layer thickness D^pc ~ 1 cited in the literature. 

The estimate of the maximum electron energy in Eq. 2.5 gives ~ 8-7 • 10*^ B^^ n_^ . 
The parameters in this expression, the magnetic field and the jet particle density, cannot be 
varied arbitrary as they also determine the jet energetics. With the anticipated model of 
the rclativistic jet dynamically dominated by cold protons, the jet (spine) kinetic power is 
given by Lj ~ it Rlj cT'j rrip (? n and the comoving magnetic field energy density should be 
less than cold protons energy density, < mpC^n. Let us rewrite the above expressions as 
Bfj^G < 2 ■ 10'^ni{^ and Lj^^j ~ 1.3 i?|^.p^r| n_4, where Lj^^j = Lj/10'^'' erg/s. Estimates often 
give for powerful jets Lj^^y ~ 1 (Sikora 2001). Hence, the number density n_4 ~ 0.01 — 1 is 
estimated for the cases of, respectively, highly relativistic and nonrelativistic large scale jets 
with Rj^kpc ~ 1 and the maximum magnetic field should be in the range 10~^ — 10~^ G. In 
Section 2 we adopted the high value of the relativistic jet density, n_4 ~ 1, as the scale of 

Q^gj g should be regarded as the lower limit for the total kinetic energy of the powerful jets 
(cf. discussion in Sikora 2001). Observations of knots in the large scale jets indicate that their 
magnetic fields amplified by the shock compression are ~ 10^^ — 10~"^ G. The magnetic field 
within the boundary layer, possibly amplified by the velocity shear, is likely to be not greater 
than these values. On the other hand, it is not expected to be smaller than the magnetic 
field ~ 10~^ — 10~^ G inside radio lobes. Therefore, we conclude that the most hkely value of 
the magnetic induction within the boundary layer is close to the assumed above B^q ~ 10. 
However, even if one will adopt the less optimistic scenario of a heavy jet (n_4 ~ 1) with the 
weak boundary layer magnetic field {B^q ~ 1) the maximum electron energy is still large, 
7eg ~ 10^, and the maximum synchrotron frequency i'syn,eq ~ -^/iG7e5 ^ 10^^ Hz. For more 
realistic parameters the synchrotron emission peaks at high X-ray frequencies. 
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As a result of weak dependence of ^eq on the boundary layer parameters the expected 
slowly varying along the jet radio-to-optical spectral index, the discussed decrease of the jet- 
counterjet radio brightness asymmetry and distributed along the jet synchrotron emission 
peaking at X-ray frequencies, are expected to hold quite generally. Presence of the flat high 
energy electron component (pile-up bump) also seems to appear in a natural way, without a 
particular tuning of the jet parameters. The main uncertainties of the presented results are 
connected with the hardly known velocity profiles of the large scale jets. The bulk Lorentz 
factor of the flow, F, determines the relative importance of the external photon fields cooling 
electrons via inverse-Compton scattering. As discussed in Section 3, the energy density of 
the AGN core emission exceeds the magnetic field energy density at distances z^pc < lO/F 
from the galactic center, while the CMB energy density is higher than ub for B^g < 3r 
in our local universe. For the considered values of 1 < F < 10 at the tens-of-kpc jet 
scales, we expect negligible contribution to the radiative output from comptonization of the 
AGN core emission, and a rough equipartition between the magnetic field and the CMB 
radiation. It results in approximately equal fiuxes of the synchrotron and the IC boundary 
layer radiation, if one neglects possible propagation effects to the observer. The form of the 
F(r) radial profile and the related spatial variation of acceleration efficiency can significantly 
infiuence the beaming pattern and intensity of the boundary layer emission. Therefore the 
large scale jet spectra at 7-ray frequencies, in particular relative normalization of the SSC 
and EC components presented at figure 2, illustrate their qualitative features only. The 
isotropic distribution of the seed photons assumed in deriving the SSC emission may result 
in its overestimate. However, even though comptonization of the CMB photons exceeds 
the self-Compton emission (for anticipated ub ^ Wcmfe ^ u^), and therefore more detailed 
treatment of the SSC radiation, involving anizotropic nature of its seed photons, does not 
result in changing of the observed composite high energy spectra. 

Finally, let us shortly comment on the jet-counterjet X-ray spectral index asymmetry 
and its angular dependence, ax,eff{(^*) (Eq- 4.2). Although the hard electron component 
is always expected to dominate the boundary layer X-ray emission, the observed spectral 
index at critical frequencies can substantially vary for different jet viewing angles, the jet bulk 
Lorentz factors and the exact shape of the electron spectrum at highest energies. Therefore, 
figure 4 illustrates only a possible behavior of the X-ray jets demonstrating different spectral 
indices of the jet and the counterjet. 
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5. Conclusions 

In the present paper, we study radiation of the uhrarelativistic electrons accelerated 
at boundary shear layers of large scale relativistic jets. The acceleration process acting at 
such layer is substantially different from the often considered first order shock acceleration. 
It can generate, in a natural way, a characteristic two component particle spectrum with a 
low energy power-law distribution extending in energy up to the harder pile-up component 
forming at higher energies due to radiation losses. In order to study a possible role of 
the mentioned process we constructed a simple, in its several aspects qualitative model of 
particle acceleration and radiation from the jet boundary. We propose, that the emission 
of such electrons can substantially contribute to the jet radiative output and we investigate 
the resulting observational effects. Modifications of the standard jet model arise from the 
kinematic effects involved (c.f. Komissarov 1990; Laing et al. 1999; Chiaberge et al. 2000; 
Celotti et al. 2001) and the mentioned nature of the acceleration process (Ostrowski 2000): 

• With the considered above continuous in situ acceleration model, the radio-to-optical 
range of the boundary layer emission is characterized by weak variation of the spectral 
index along the jet, assuming the boundary layer physical parameters change smoothly 
with the distance from the nucleus. The considered radial velocity profile of the rel- 
ativistic jet defines the beaming pattern of the boundary layer emission. In general, 
the velocity shear within the boundary region results in decreasing the observed jet- 
counterjet brightness asymmetry, affecting evaluations of the jet bulk Lorentz factors. 

• If the considered cosmic ray population involves the high energy bump electrons with 
the energy density comparable to the magnetic field energy density, one should observe 
a high frequency (X-ray) excess over the power-law radiation extrapolated from the 
lower frequency measurements. In a framework of our model, the spectral properties 
of the X-ray bump emission (the effective spectral index, jet-counterjet brightness 
asymmetry, etc) can substantially differ from the one observed at lower frequencies. 
We speculate, that such features can be seen at least in some Chandra jets. 

• The considered model can explain the observed in some cases edge-brightened jets in 
a natural way. 

One should note, that our conclusions depend weakly on the unknown detailed physical con- 
ditions within the transition layer. The model presented here, far from being complete, is 
intended to describe the main radiative properties of the boundary region under simple as- 
sumptions about turbulent character of the shear layer, relativistic jet velocities on the large 
scales and energy equipartition between the magnetic field and the accelerated electrons. We 
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expect, that careful analysis of observational data should reveal the effects discussed above 
and thus provide restrictions both on physical parameters within such boundary layers and 
the proposed model of electron acceleration. 

We are grateful to Marek Sikora for his help and discussions. MO acknowledges a useful 
discussions with Vahe Petrosian. Remarks of the anonymous referee enabled to substantially 
improve the final version of the paper. The present work was supported by Komitet Badah 
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Fig. 1. — A simple model spectrum A'"(7) = d ne{'y)/d log 7 of particles accelerated within 
the turbulent shear layer. The acceleration process creates power-law energy distribution 
ended by the cut-off (the dashed line) or by the high energy pile-up bump at the maximum 
electron energy (the soild hne) (cf. Ostrowski 2000). 
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Fig. 2. — The observed spectral energy distributions of the radiation generated within the 
boundary shear layer for = 10^ and 9* = 5^ , 30° , 60° , 90°, at the respective panels. The 
presented spectra are integrated over the flow Lorentz factor profile (3.9) and correspond to 
the parameters discussed in the text. The index '1' refers to the emission oc of the 

power-law part of the electron spectrum. The index '2' refer to the emission of the high 
energy electron bump component. In the case of the SSC process, monoenergetic electrons 
contribute also to the emission oc Absorption of VHE 7-rays during the propagation 

to the observer is neglected. 



-29- 




Fig. 3. — A ratio of the observed radio fluxes from the jet and the counterjet, S*{9*)/S*{n — 
^*), as a function of the viewing angle 9*. The dotted hues correspond to the model with the 
uniform flow Lorentz factor Tj = 10 or 3 within the whole jet. The solid lines correspond to 
the radiation from a boundary shear layer with the radial profile (3.9) and the same Fj — 10 
or 3. The Lorentz factors are provided near the respective curves. 
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Fig. 4. — An illustration of the possible effective X-ray spectral index of a jet ('j') and a 
counterjet ('cj') boundary layer emission, ax,eff, as a function of the viewing angle 9*. The 
solid lines correspond to Tj — 10 and the dashed ones to Tj — 3. 



